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ABSTRACT 

We study torsional Alfven oscillations of magnetars, i.e., neutron stars with a strong 
magnetic field. We consider the poloidal and toroidal components of the magnetic 
field and a wide range of equilibrium stellar models. We use a new coordinate system 
{X,Y), where X = ^/oisin^, Y = ^^cos^ and ai is the radial component of the 
magnetic field. In this coordinate system, the 1 + 2-dimensional evolution equation 
describing the quasi-periodic oscillations, QPOs, sec Sotani ct al. (2003), is reduced 
to a 1 + 1-dimensional equation, where the perturbations propagate only along the 
y-axis. We solve the 1 + 1-dimensional equation for different boundary conditions and 
open magnetic field lines, i.e., magnetic field lines that reach the surface and there 
match up with the exterior dipole magnetic field, as well as closed magnetic lines, 
i.e., magnetic lines that never reach the stellar surface. For the open field lines, we 
find two families of QPOs frequencies; a family of "lower" QPOs frequencies which is 
located near the X-axis an d a family of "upper" frequencies located near the F-axis. 
According to iLevinI (j2007[ ). the fundamental frequencies of these two families can be 
interpreted as the turning points of a continuous spectrum. We find that the upper 
frequencies are constant multiples of the lower frequencies with a constant equaling 
2n + 1. For the closed lines, the corresponding factor is n -t- 1 . By these relations, we 
can explain both the lower and the higher observed frequencies in SGR 1806-20 and 
SGR 1900-f 14. 
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1 INTRODUCTION 

Strongly magnetized compact stars, so called 'magnetars', 
reveal their presence via giant flares with peak luminosities 
of 10*'* — 10*® ergs/s. The giant flares are typically accom- 
panied by a decaying tail which may last several hundreds 
of seconds. These sources are known as Soft Gamma Re- 
peaters (SGRs) , and up to now, three giant flares have been 
detected, SGR 0526-66 in 1979, SGR 1900-hl4 in 1998, and 
SGR 1806-20 in 2004. The analysis of the timing of the 
last two events revealed several QPOs in the decaying tail 
whose frequencies are approximately 18, 26, 30, 92, 150, 625, 
and 1840 Hz for SGR 1806-20, and 28, 53, 84, an d 155 Hz 
for SGR 1900-hl4, see I Watts fc Strohmaved (|2006D . It is be- 
lieved that during an SGR event torsional oscillations in the 
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solid crust of the star could be excited, (|Duncanlll99i ). that 
lead to the observed frequencies in the X-ray tail. Indeed, the 
frequency of many of these oscillations fit the values of the 
torsional mode oscillations of the solid crust of a compact 
star. However, since not all of the obs erved frequen c ies ar e 
explainable by pure crust oscill ations, ISotani et al.l l|2007h : 
ISamuelsson fc AnderssonI (|2007l ), alternative scenarios have 
also been suggested. For example. iGlampedakis et al.l (j2006l ) 
claimed that the observ ed spectra m ay be explained via 
Alfven oscillations, while iLevinI (|2006l ) pointed out that the 
Alfven oscillations of magnetars co uld be conti nuous. Re- 
cently, a toy model calculation by iLevinI (|2007l ) indicated 
that the edges or turning points of the continuum could 
correspo nd to long-lived QPO s, and according to Levin, the 
results of lSotani et al.l (|2008al ) indicated for realistic magne- 
tar configurations that the Alfven oscillations of magnetars 
could be continuous which would explai n the lower observed 
frequencies. In a more recent paper by ISotani fc KokkotasI 
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l|2009l ). it was shown that the spectrum of the polar Alfven 
oscillat ions is actual l y discre te. 

In ISotani et al ] (l2008a'), the authors performed two- 
dimensional numerical simulations of linearized Alfven oscil- 
lations in magnetars. Their model improves the previously 
considered toy models in various ways. General Relativity 
is assumed, various realistic equations of state (EOS) are 
considered and a consistent dipolar magnetic field is con- 
structed. However, it does not take into account the pres- 
ence of a solid crust and only examines the response of the 
ideal magnetofluid to a chosen initial perturbation. The two- 
dimensional partial differential equation (PDE) that was 
used to study the Alfven oscillations has a mathematical 
pathological behavior, as will explained later, that forced 
them to introduce artificial numerical viscosity in an at- 
tempt to stabilize the numerical evolution. The numerical re- 
sults presented in Sot ani et al.l ([2008a) are compatible with 
the observations. For example, they found two families of 
QPOs, corresponding to the edges or turning points of a 
continuum, with harmonics at near integer multiples. With 
this identification, they could set an upper limit to the dipole 
magnetic field of ~ 3 to 7 x lO^^G, and they could limit 
the models to very stiff EOS for values of the magnetic field 
strength near to its upper limit, or moderately stiff for lower 
values of the magnetic field. 

In an extension of that paper, ISotani et all (|2008bh . 
studied the axisymmetric crust torsional modes of magne- 
tars with poloidal and toroidal magnetic field, where both 
components are confined in the crust. The numerical results 
showed that this magnetic configuration cannot explain the 
actual observational da ta of SGRs, and is in agreem ent with 
a more recent result by Ivan Hoven fc LevinI l|2008l ). 

After a more elaborate study of the master equation, 
we discovered the origin of the numerical instability, and, 
moreover, we managed to reduce it, by using a coordinate 
transformation, to a 1 -I- 1-dimensional PDE which can be 
evolved stably without the need of artificial viscosity and 
also allows a semi-analytically estimate of the properties of 
the continuum and the associated QPO's. 

In this paper, we adopt the units of c = G = 1, where c 
and G denote the speed of light and the gravitational con- 
stant, respectively, and the metric signature is (— , -f ). 



2 THE EQUATION OF THE PROBLEM 

We consider a spherically symmetric and static star de- 
scribed by the TOV equations and the line element 



dt + e 



2A(r) 



dr^ + r^{de^ + sm{e)d(t)'^) . (1) 



The structure of the star is not infiuenced by the presence of 
the magnetic field because the last has an energy Em which 
is a few orders of magnitude smaller than the gravitational 
binding energy Eg. Typically, E^/Eg ~ 10-*(B/10"G)^ 

Generally, spherical stars show two type of oscillations, 
spheroidal oscillations with polar parity and toroidal oscil- 
lations with axial parity. The observed QPOs in SGR X- 
ray tails may originate from toroidal oscillations, since these 
can be excited more easily than poloidal oscillations because 
they do not involve density variations. 

The MHD oscillations of the equilibrium model are de- 
scribed by the linearized equation of motion and by the mag- 



netic induction equations, see ISotani et all |20o3). These 
two perturbation equations can be converted into a, 1 + 2- 
dimensional evolution equation for the displacement func- 
tion y{t,r,0) which is related to the contravariant compo- 
nent of the perturbed 4-velocity, Su''' via 

Su'^ = e-'^dtyit,r,e). (2) 

For an analytical derivation, see ^Sot ani et al ] (|2008a'). The 
1 + 2-dimensional time-evolution equ ation for the displace- 
ment y{t, r, 9) has the following form l|Sotani et al.l (|2008al )): 



A, 



. d^y , d^y . d^y . dy , dy 

(3) 

All the coefficients Au, A20, A\\, Ao2-, A\o and Ao\ depend 
on the coordinates r and 9, but not on time. They are given 
by 

^2 
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where ai(r) is the radial component of the electromagnetic 
four-potential and ji{r) is the four-current. A prime indi- 
cates a derivative with respect the radius. The symbol /i 
denotes the shear modulus. In the following, we will set ^ 
equal to zero, neglecting the presence of a solid crust, since 
we will focus on the oscillations of the magneto-fiuid core 
that plays a m ajor role in m agne tars QPOs. In fact, as it ha s 
been shown bv lLevinI (120061 ) and lGlampedakis et ahl (|2006l ). 
in the presence of a strong magnetic field pure crustal modes 
cannot explain the present QPOs frequencies. As a conse- 
quence, it is necessary to study global modes, i.e., Alfven 
waves propagating through the whole star. 

The radial function ai is determined by solving the 
Grad-Shavranov equation: 

2 



ai e 



+ ($' 



a/\ -2A / 

Ale ai 



-47rji 



(10) 



As boundary conditions, we assume regularity at the center, 
ai ~ aor^, where ao is a constant, and continuity at the 
surface with a di pole field outside th e star. 

According to lSotani et al.l (|2007l 'l. the solutions of equa- 
tion (131) have to satisfy the following boundary conditions: 

• regularity at the center: ^ = at r = 0; 

• no traction on the surface: y^r = at r — R; 

• axisymmetry: y^g — at 9 = 0; 

• equatorial plane symmetry for £ — 2 initial data: 3^ = 
at 9 = n/2; 

• for £ — 3 initial data: y,e — at 9 — -k 12. 
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Figure 1. Plot of a\(r) for different EOS. The coordinate r is 
normalized to the stellar radius. 



In lSotani et al.l l|2008al ). in order to avoid numerical insta- 
bility, equation ((Sj was solved by introducing an artificial 
4-th order 'dissipation' term. Here, instead, we find a coor- 
dinate transformation that allows us to avoid the use of a 
dissipation term and, as we will show later, to avoid the sin- 
gular behavior of the time independent form of the equation 
((3|. In fact, if, in a first approximation, we neglect the first 
derivatives in the direction of r and 9 in the equation @, 
then only the coefficients A2Q, A\i and A02 remain. Calcu- 
lating the determinant of the principal part of the equation 
(121, we find that 



D = A20A0 



2 12 



An_ 
2 



cos^(6l)sin^(6') 



2 12 



cos^(e')sin^(e») = 



(11) 



Hence the equation is parabolic in every point of its domain, 
and therefore does not describe propagation of 2D waves. 
Significant simplification is achieved by introduction of a 
new coordinate system given by the transformation 



X = 



/fli sm f 



Y = 



/fli cos t 



(12) 



This new coordinates system is a kind of Cartesian coordi- 
nate system where, instead of the radial coordinate r, we use 
the radial component of the magnetic four-potential ai . The 
behavior of ai inside the star is determined via the equa- 
tion (|10p . On the surface, we require the continuity of ai 
with a dipole solution outside. In Figure[TJ we plot the func- 
tion ai(r) for different equations of state (EOS). The sub- 
scripts indicate the mass of the star. For example, WFF14 
corresponds to an equ ilibrium stellar model with EOS WFF 
l|Wiringa et al.l (|l988l )) and M = IAMq. 

The magnetic field and the radial function ai are 
related by 

= e^f'^-'u^F^, , (13) 

where e"^'^^ is the Levi-Civita tensor, up is the background 
four-velocity and is the electromagnetic tensor. Then 
equation (|13|l leads to 



Ur = _ cos f 



Ha = 



4-K 



(14) 



where we define Ha ~ Bajv^- Note, however, that the co- 
ordinates defined by (|12p are not the usual ma gnetic coordi- 
nates known in the literature, see for example I Alladio et al.l 
( 199 6) for a review, since they mix the poloidal angle and 
the radial component of the magnetic field in (X, y). 

When we transform equation (|3} according to (|12p . we 

get: 
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(15) 

By substituting a\" from equation (|10|) into the second and 
third term on the right hand side of equation (jlSp , we arrive 
at 

d^y 



Au 



T d^y,2 dy 



where 

Aio -- 

Aoi -- 
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2-Kr- 



■ai 



2 2A 



A ■ 2 
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(16) 

(17) 
(18) 



i.e., the coefficient of dy/dX vanishes identically. The equa- 
tion p6|) is a 1 + 1-dimensional wave equation. The pertur- 
bations propagate only along the Y direction. In the last 
direction, the problem reminds of a study of waves propa- 
gating along the strings of a musical instrument. 

The construction of the grid in these {X, Y) coordinates 
needs special attention. The function ai has a maximum in- 
side the star, and the position of this maximum depends on 
the strength and the topology of the magnetic field. A strong 
magnetic field, and also a toroidal field added to the stan- 
dard poloidal field, pushes the maximum near the surface of 
the star. The position of ai max determines also the number 
of the closed magnetic field lines confined inside the star. 
Since the closed magnetic field lines don't reach the surface, 
they are subject to different boundary conditions than the 
open magnetic field lines. This difference leads to a different 
behaviour of the oscillations of open and closed magnetic 
lines, as will be shown later. In practice, we split the prob- 
lem into two parts for both cases, i.e., for closed and open 
magnetic field lines. For the open magnetic field lines, the 
perturbation function u{X, Y) is evolved until some points 
Xmax and Ymax that correspond to the maximum of ai max, 
i.e., we work on the part A of Figure [2l Then the evolution 
continues, in the part Ai of Figure [2l from the points Xmax 
and Ymax and ends at the points Xsurt and Vsurf that cor- 
respond to the surface of the star. Note that in this second 
part, we have to use the transformation 



X^-^smO y = -Va7cos6l, (19) 

due to the sign change of the derivative of oi. In a similar 
fashion for the closed magnetic lines, we evolve u{X, Y) until 



Xmax and Ym 



but then the simulation doesn't end on the 



surface, but at Xeq and Yeq, that correspond to points on 
the equator. From a technical point of view, we just '^operi' 
the field lines, then we solve the two problems separately 



4 A. Colaiuda, H. Beyer & K. D. Kokkotas 




Figure 2. Plot of the magnetic field lines in the (X,Y) coordi- 
nates. Note that, after the point X{ai max) and Y{ai max): lines 
of the same color continue in different quadrants. For example, 
the lines of the domain A continue in the domain Ai. 



and match the two solutions that we get via the boundary 
conditions. This procedure will be explained later. Figure [2] 
may help in understanding the procedure followed in solv- 
ing this problem. The equation (|16|) must be solved with the 
appropriate boundary conditions. The boundary condition 
are different for open and closed magnetic lines. The open 
lines reach the surface and match up with the dipole mag- 
netic field outside the star. Then, as we will show later, they 
couple on the surface and oscillate together. Differently, the 
closed magnetic lines confine themselves in the interior of 
the star, forming closed loops, and they never reach the sur- 
face. So they don't have any coupling to the open magnetic 
lines and oscillate independently. When we 'open' the closed 
magnetic lines, creating magnetic strings, we open the closed 
loop, but don't change the boundary conditions at the ends 
of the strings. 

In the {X, Y) coordinate system, the boundary condi- 
tions become 

• regularity at the center: y{X, Y) = for X = and 

y = 0, 

• no traction on the surface for the open lines: 



a-L 
2ai 



dy dy 



= 0; 



(20) 



axisymmetry at X = 0: 



Y 



dy 

dX 



0; 



• equatorial plane symmetry for 

y{x,Y) = at y = 0. 

• for £ = 3 initial data at y = 0: 



y 



dy 
dx 



0; 



(21) 

= 2 initial data: 
(22) 



As a consequence, for initial data with 1 = 2, the open 
'magnetic strings' are coupled on the surface of the star and 
near the axis for X = Q (corresponding to 9 = 0). With 
initial data 1 = 2), there is an additional coupling at the 
equator (corresponding to y = 0). 



2.1 Numerical method 

We study a variety of equilibrium models, constructed for 
three representative EOSs and different masses that range 
from 1.4 — 2.0 Mq, see Table [1] We set the magnetic field 
strength to = 4 x lO'^^ Gauss, matching the field strength 
extrapolated from the analysis of the observations for SGR 
1806 - 20 and SGR 1900 -I- 14. We construct a numerical 
equidistant grid 90 x 90 in the [X, Y) coordinates, setting 
Xmax = V^i max siu 6^ and y^ax = max COS 6 and Vary- 
ing 6 from and 7r/2. The structure of our grid is displayed 
in the first quadrant of Figure [2] Note that each 'magnetic 
string' in the first quadrant continues in the third quadrant. 
We evolve these two parts in time like two separate prob- 
lems, taking into account the different boundary conditions 
for open and closed magnetic field lines. The two parts are 
linked in every time step, requiring the continuity of the 
function 3^ and its derivative 3^' . For our evolution, we use 
an iterative Crank-Nicholson scheme and we evolve the per- 
turbations for more than Is (see Figure In order to check 
the stability of the scheme, we performed an evolution for 
2s and found that our scheme was still stable. 

We evolve the perturbations by providing initial data 
corresponding to either odd, ^ = 3 or even, 1 = 2, boundary 
conditions at y = 0. A change in the value of I leads to 
different boundary conditions for open and closed lines and 
leads to different scaling in the oscillation frequencies of the 
open field lines. 

We compute the FFT of various 'strings' and find 
that each 'string' c orresponds to a different frequency. In 
ISotani et al.l l|2008al ). the authors find a family of upper fre- 
quencies near the magnetic axis for 9 = Q, corresponding to 
points along the Y-axis in our coordinates, and a family of 
lower frequencies for 9 = 12, corresponding to points along 
the X-axis in our coordinates. They identify the fundamen- 
tal frequency peaks of the two families of QPOs with the 
ed ges or t urnin g points of the Alfven continuum suggested 
bv lLevinI (|200'if ). These two frequencies were named funda- 
mental lower Lo, for 9 = 7r/2, and fundamental upper C/q, 
for 9 = Q QPOs. The overtones are indicated by L„ and i7„. 
In the following, we adopt these conventions when we refer 
to the oscillations of the open field lines, while we call the 
frequencies of the closed field lines C„. 

Ou r results confirm partially the results of ISotani et al.l 
(l2008al ). For open lines, a family of upper frequencies ap- 
pears near the Y-axis, while a family of lower frequencies 
appears near the last open magnetic field line, just before the 
region of closed lines. This can be explained by the fact that 
the 'strings' near the Y-axis correspond to a stronger average 
magnetic field because in the (r, 9) coordinates they corre- 
spond to denser field lines, see Figure [G] Since the frequency 
of the Alfven oscillations is proportional to the strength of 
the magnetic field, we expect that these frequencies are the 
upper ones. In addition, we find that, for all the 'magnetic 
strings', the upper frequencies are multiples of the funda- 
mental frequency both in the case of I/„ and in the case of 
[/„, see Figure |4] However, this multiplicity is different when 
we consider even t = 2 ax odd £ = 3, initial data. In fact, we 
obtain 



(2n + 1)/lo 
(2n + l)/c„ 



(23) 
(24) 
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/^odd ~ (^^ + 1)^0 (25) 
/^odd ~ {n + l)fu,, (26) 

where the superscripts even and odd indicate resuhs derived 
for even 1 = 2 and odd ^ = 3 parity initial data and bound- 
ary conditions on the equator. The different scaUng shown 
in the relations ((23]) - (HU and ([251)- ((26l) is due to the dif- 
ferent boundary conditions used for odd and even values of 
I. 

The problem can be understood by looking at a simple 
problem of two strings of finite length connected by by an- 
other spring. If the two strings are connected at just one of 
their ends, while the others are kept fixed, like in the case 
of even parity where the 'magnetic strings' are coupled just 
near the surface and are fixed at the equator, then we obtain 
the relation s (12311 and (I24II. Thi s type of configuration was 
analyzed in iKokkotas fc Schut j (Il986l ). in order to explain 
the w-modes of neutron stars. An analytic derivation of the 
relations (|23p and (|24|) can be found there. If instead the 
two strings are coupled in both ends, like in the case of odd 
parity, then the relations (|25p - (|26|) follow. 

The frequencies of the closed magnetic field lines show 
a different behaviour. The upper and lower frequencies are 
scaling as, see also Figure O 

/c„ ^ (n + l)/co . (27) 

The difference between the relation (|23|l and (|27p is due 
to the different boundary conditions. The closed lines never 
reach the surface. They start from the equator where they 
close. This means that they are subject to the same bound- 
ary condition at both of their ends which resembles the toy 
problem of a string that oscillates with both ends fixed and 
gives the n-\-l behaviour, as in the case odd values of I for the 
open field lines. The frequencies that we find for closed lines 
initially decrease to a minimum and from then on increase. 
This behaviour results from two factors. First, from the fact 
that closed lines are shorter than open lines, and, second, 
from the fact that the magnetic field becomes somehow lo- 
cally weaker as we move off the axis. This behaviour of the 
spectrum has been identified after comparing our data with 
th e results coming from the non-linear evolutions presented 
in lCerda-Duran et all (|2009l ). 

The relations here provide the limits of the continuous 
part o f the spectrum, in accordance with the results of lLevinI 
l|2007l ). Our results are summarized in Table [T] 

2.2 Why we don't need artificial viscosity 

In the preceding work lSotani et al] (|2008al ). an artificial vis- 
cosity was added in order to avoid a fast growing instability. 
This instability is du e to the presence of the so called " O- 
type neutral point" l|Tavloij (| 19791 ) V If the magnetic field 
lines of a poloidal magnetic field are closed and confined to 
the stellar interior, then there is a point on which these mag- 
netic field lines converge and become degenerate (see Figure 
[G]). This point corresponds to the maximum of the function 
ai. 

In our scheme, the magnetic field lines are practically 
straight, see Figure [2j even in the vicinity of the point that 
corresponds to ai max. This behavior is due to our choice to 
'open' the grid and to evolve the open and closed field lines 
as two independent problems. When we look at the point 
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Figure 3. The time evolution of dty for an open line in a point 
near the Y-axis for the stellar equilibrium model WFF14 and 

C = o. 
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Figure 4. The FFT of a point near the Y-axis, green/dashed line, 
and near the critical point, red/continuous line, for open lines and 
even parity, £ = 2. The stellar model is WFF14. 
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Figure 5. The FFT of the first closed line for the stellar equilib- 
rium model WFF14. It is worth noticing the spacing of the peaks 
in accordance to I I27II 
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Figure 6. The magnetic field fines for the equilibrium steUar 
model WFF14 and C = 0. 



0.1 max, we find that it practically doesn't evolve and, because 
the 'strings' are coupled only along the surface and near the 
X = axis, cannot be influenced by the oscillations of the 
other 'strings'. In this way, we avoid the singular behaviour 
shown in the (r, 6) coordinates. The last is the main reason 
why our code is stable also without artificial viscosity. 



2.3 Toroidal fields 



From equation (|13|l , the toroidal magnetic field can be writ- 
ten as 



An 



sin 



(28) 



where is a constant only depending on the toroidal mag- 
netic field. By increasing or decreasing we can infiuence 
the strength of the toroidal magnetic field. So ( can be 
interpreted as a parameter that describes the strength of 
the toroidal magne ti c field in relation to the poloidal field, 
see IColaiuda et al.l l|2008l ) for details. The presence of a 
toroidal field modifies equation (1101) by adding an extra term 
'^"ai (see lSotani et al.l l|2008bt )). 



ai e 

(29) 

Hence, the toroidal field influences our evolution equation 
(|16|) only through the presence of the extra term ai" . By 
solving equation (|29|l for ai" and subsequent substitution 
into equation p6p . the toroidal magnetic field appears in 
the 1-1-1 equation (|16|) only through the parameter in the 
term 



7 X 
Aoi = yai 



2e2A„j 



2 -2# 0,1 . 

C oie 4TVJ-, 



(30) 



This means that the toroidal field affects the local propa- 
gation speed of the Alfven modes. Also, the toroidal field 
pushes the position of ai max outwards and contributes to 
the stability of the magnetic field in the star. In Figure [T] 
we plot the function ai for the EOS WFF14 and for different 
values of the parameter It is clear that when ( increases, 
i.e., the strength of the toroidal magnetic field increases, 
then ai max is pushed closer to the surface . 

Hence, it is natural to expect changes in the oscillation 



WFF,4C=0 - 
WFF14 ^=0.18 
WFFia C=0.20 - 
WFF14 i;=0.24 



Figure 7. The function ai as a function of the normalized coordi- 
nates r/R for an equifibrium model constructed with the equation 
of state WFF and a mass M = IAMq. Different curves reflect 
different values of the parameter i.e., different strength of the 
toroidal magnetic field. 



spectrum due to the presence of the toroidal component of 
the magnetic field. We find that the frequencies are lower 
compared to the case in which only a poloidal field is present, 
see TableO especially for the closed magnetic 'strings'. This 
can be explained because, by the presence of a toroidal mag- 
netic field ai max, i.e., the point were the magnetic field lines 
converge, moves outwards, and hence the 'magnetic strings' 
increase in length. Since the frequency is inversely propor- 
tional to the length of the 'strings', as the length increases, 
the frequency decreases. 



3 IDENTIFICATION OF THE QPOS 
FREQUENCIES 

In a paper bv I Watts fc Strohmaveij l|2006l ). the timing anal- 
ysis of the SGR 1806-20 and the SGR 1900-M4 has been 
computed. The study shows the appearance of several QPOs 
in the tail of these events. In particular, for SGR 1806-20, 
the identified frequencies are 18, 26, 30, 92, 150, 625 and 1840 
Hz, while for SGR 1900+14, the frequencies are 28, 53, 84 
and 155 Hz. 

Using the relations (|25[) and (1261) . we can easily show 
that for SGR 1900+14, if we consider /^odd = 28Hz, we can 
identify /j^odd = 56Hz, and /^odd — 84Hz. In a similar way 
for SGR 1806-20, we have ^gvd = 18Hz and fjjadd = 30Hz, 
and consequently /[7|von ~ 92Hz or fjjodd — 92 , fuf"^ — 
150Hz or /[/odd ~ 150. If we look at our results in the Tables 
1 and 2, we can identify the QPOs in SGR 1900+14 with the 
lower frequencies for odd parity of a model with equation 
of state WFF and a mass around I.8M0. For the QPOs 
from SGR 1806-20, we can fit the observational data by 
a model with equation of state APR, mass around 1.4Mq 
and a toroidal field a little stronger than the poloidal field , 
C = 0.18. 
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Table 1. Frequencies of Alfven QPOs for a representative sample of equilibrium models, constructed for different EOS and fixed magnetic 
field strength S = 4 X 10^^ Gauss. The frequencies are computed for points near the critical point, L„, and near the Y-axis, Un, for 
open lines, both for even and odd parity, and for the closed line (Cn)- 
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69.55 


93.11 


84.15 
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116.1 
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0.264 
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54.41 
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66.75 
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67.31 


81.89 


80.21 


99.84 


42.06 
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47.06 
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24.04 
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95.17 


86.28 


121.8 


74.78 
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142.2 
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180.0 









19.04 


39.68 


21.69 


44.97 


19.57 
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58.72 


79.88 


66.6 


89.40 


39.68 
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96.28 


119.6 


109.5 


133.8 


60.31 
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15.25 


31.60 


13.07 


APR20 
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37.59 


50.09 


44.12 


61.56 


25.60 
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62.10 


74.63 


74.63 


92.61 


38.13 









13.10 


27.83 


16.78 


35.61 


13.52 




0.198 


1 


40.11 


55.66 


51.16 


70.40 


27.85 
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68.35 


83.50 


85.54 


105.6 


41.34 



Table 2. Frequencies of Alfven QPOs and their ratios for different stellar models with a toroidal component of the magnetic field and a 
surface magnetic field _B = 4 X 10^^ Gauss. For open magnetic 'strings' the frequencies are taken for points near the critical point, Ln, 
and near the Y-axis, Un, both for even and odd parity. We also give the value of the frequencies for the lower closed magnetic 'string' 
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..{Hz) 
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15.15 


30.30 


16.24 


33.94 


15.76 


WFFn 


0.24 
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46.06 


61.20 


50.91 


66.66 


31.15 
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75.75 


90.90 


81.20 


99.99 


79.38 









18.95 


38.46 


22.85 


49.6 


19.51 




0.18 


1 


56.83 


76.91 


69.67 


98.65 


39.01 
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95.30 


114.8 


117.0 


147.7 


57.96 









16.21 


32.42 


19.35 


41.84 


16.73 


APRu 


0.20 


1 


48.11 


64.84 


59.05 


83.67 


33.47 
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64.84 


97.79 


99.88 


125.5 


51.25 









11.47 


22.93 


13.92 


29.90 


11.88 




0.20 


1 


34.40 


45.46 


42.18 


58.57 


23.75 






2 


56.52 


68.39 


70.05 


88.46 


35.63 



4 DISCUSSION AND CONCLUSION 

We have investigated torsional Alfven modes of relativis- 
tic stars with poloidal and toroidal magnetic field using a 
new coordinate system {X,Y), where X = -^aT sin 6* and 
Y — ,/aicos9. In this new coordinate system, the per- 
turbed two dime nsional equat i on for the toroidal displace- 
ment y{t, r, 9) in lSotani et al.l (|2008al ) is reduced to a 1 + 1- 
dimensional equation, and the perturbations, i.e., Alfven 



waves, are propagating only along the Y-direction. By solv- 
ing the problem with the appropriate boundary conditions 
for open and closed magnetic 'strings', we find that each 
'magnetic string' corresponds to a family of frequencies. The 
upper frequency belongs to the family of open magnetic 
'strings' near the Y-axis because the Alfven frequency is 
proportional to the local strength of the magnetic field, and 
since near the Y-axis, the magnetic field is stronger than in 
the rest of the star. The family corresponding to the lower 
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frequency is the one just before the 0-point, i.e., the point 
where the magnetic Unes converge and degenerate into a 
point. We find that for each family of frequencies there is 
a relation between the overtone frequencies and the funda- 
mental frequencies, and this relation changes if we consider 
odd or even parity, see equations p3|l . (f24|l . ([25}, p6|l . In 
addition, we find a third family of frequencies linked to the 
closed lines. This family shows a minimum, whose value is 
similar to that of the lower frequencies of the open lines 
for even parity. Also, we find that in presence of a toroidal 
magnetic field the frequencies are lower than in the case 
with only a dipole magnetic field, both for closed and open 
magnetic 'string'. This is due to the change in the position 
of oi max that moves outwards. This results in longer 'mag- 
netic strings' and, as consequence, to a lower characteristic 
frequency for each 'string'. 

In a future work, it will be interesting to study whether 
the presence of a solid crust can quantitatively and qualita- 
tively alter our results. 
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